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We report a gate-controlled transition of a semimetallic InAs/GaSb heterostructure to a topolog¬ 
ical insulator. The transition is induced by decreasing the degree of band inversion with front and 
back gate voltages. Temperature dependence of the longitudinal resistance peak shows the energy 
gap opening in the bulk region with increasing gate electric field. The suppression of bulk conduc¬ 
tion and the transition to a topological insulator are confirmed by nonlocal resistance measurements 
using a dual lock-in technique, which allows us to rigorously compare the voltage distribution in the 
sample for different current paths without the influence of time-dependent resistance fluctuations. 


Topological insulators (TIs) have attracted strong in¬ 
terest as a new quantum state of matter categorized nei¬ 
ther as a metal nor an insulator.—The unique elec¬ 
tronic properties of TIs originate from their band struc¬ 
tures characterized by conduction-valance band inver¬ 
sion with an energy gap opening in the bulk region. 
The band inversion results from the interplay between 
the alignment of bands with different orbital characters 
and the strong spin-orbit interaction inherent to materi¬ 
als containing heavy elements, such as Bi compounds^— 
and HgTe.— - — Owing to their topological nature, gap¬ 
less states with distinct transport characteristics emerge 
at the surfaces and edges of three- and two-dimensional 
(2D) TIs, respectively. In a 2D TI, also known as a quan¬ 
tum spin Hall insulator, the edge state comprises coun- 
terpropagating channels with opposite spin— Since 
elastic back scattering is prohibited by time-reversal 
symmetry, dissipationless spin-polarized transport is ex¬ 
pected. The observation of the quantum spin Hall effect 
in a HgTe/CdTe quantum wells - — has established it as 
a prototypical system for 2D TIs. 

Recently, InAs/GaSb heterostructures, with the so- 
called “broken-gap” or “type-II” band alignment, have 
attracted interest as a route to achieving a 2D TI by 
a combination of materials each having a trivial band 
structure typical of III-V semiconductors— - — Different 
from conventional TIs whose topological properties orig¬ 
inate from the band structure of individual materials, 
the band inversion in InAs/GaSb is provided by the rel¬ 
ative band alignment of the two semiconductors, where 
the InAs conduction-band bottom is located below the 
GaSb valence-band top. Experiments have so far con¬ 
firmed the existence of edge channels in InAs/GaSb het- 
erostructures in the band-inverted regime— - — Addition¬ 
ally, it has been shown in Ref. [U that the band alignment 
can be tuned by the thickness w of the InAs layer. This 
is manifested by the transport behavior’s changing from 
that of a normal band insulator through that of a TI to 
semimetallic as w is varied from 10 to 14 nm. However, 
the major advantage of the InAs/GaSb system -the band 
alignment can be controlled in situ by gate voltage s 11 ’ 17 
has not yet been fully exploited. 


In this paper, we study the effects of gate electric field 
on the transport properties of micrometer-sized dual-gate 
InAs/GaSb heterostructure device using nonlocal trans¬ 
port measurements. Different from HgTe/CdTe, in which 
both the electron and hole wave functions are confined to 
HgTe and the band inversion is determined by the HgTe 
layer thickness, the electron and hole wave functions in 
InAs/GaSb are confined to different layers, which allows 
us to tune the degree of band inversion in situ by ap¬ 
plied gate voltages—We demonstrate that an origi¬ 
nally semimetallic heterostructure can be tuned into a 
TI using gate voltages. Our results represent an impor¬ 
tant aspect of the InAs/GaSb system, which is useful for 
studies of TIs. 

Figure 1(a) shows schematically the band diagram of 
the system we study in this paper. The InAs/GaSb het¬ 
erostructure is sandwiched between Al^Gai-^Sb barriers 
and has gate electrodes on both sides. We focus on the 
band-inverted regime, that is, the situation where the 
thicknesses of the InAs and GaSb layers are such that 
the bottom of the electron subband in InAs (with energy 
E e o) is located below the top of the hole subband in GaSb 
(with energy Eho) as shown in Fig. 1(c). The electron and 
hole wave functions, which are predominantly confined 
to the InAs and GaSb layers, respectively, are hybridized 
through the heterointerface to produce an anticrossing 
gap A at a finite in-plane wave vector k [Fig. 1(b)]. As 
a consequence of the inverted band alignment, linearly 
dispersing bands connecting the upper and lower bands 
emerge [dotted lines in Fig. 1(b)], which correspond to 
the edge channels. When the Fermi level is in the energy 
gap, the system becomes a 2D TI, where the conductivity 
in the bulk region is expected to vanish and the transport 
is governed by the edge channels. Experimentally, while 
the existence of the anticrossing gap has been confirmed 
by capacitance spectroscopy^ transport measurements 
usually indicate a resistance value that remains fairly 
low when the Fermi level is swept across the expected 
gap region— - — One possible origin of this behavior is 
the semimetallic band structure that results from the 
anisotropy of the GaSb valence band^ [Fig. 1(e)], which 
we will discuss later. 
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FIG. 1. (Color online) Schematics of (a) band-edge profile of 
the sample used and (b) dispersion relations in a TI. (c) and 
(d) Schematics illustrating the effects of gate electric field on 
the squared wavefunctions of electrons (red) and holes (blue) 
(upper panels) and band-edge profile (lower panels). The di¬ 
agrams in (c) and (d) represent the flat-band condition and 
the case with the gate electric field increasing the electron- 
hole wavefunction overlap, respectively. E e o and Fho indicate 
the energies of the electron and hole ground subbands, re¬ 
spectively. (e) and (f) Schematic illustrations of effects of 
gate electric field on the dispersion relation in the presence of 
valence-band anisotropy, (e) and (f) correspond to the situa¬ 
tion in (c) and (d), respectively. The shaded regions indicate 
the energy range of the anticrossing gap for each k direction. 


Figures 1(d) and (f) illustrate the effects of an electric 
field on the band structure. Applying a positive voltage 
on the GaSb side decreases E^q while a negative voltage 
on the InAs side increases E e o- This results in a lesser 
degree of band inversion as measured by E^q — E e o. Con¬ 
sequently, the anticrossing is shifted to a smaller k and, 
at the same time, the larger spatial overlap of the elec¬ 
tron and hole wave functions leads to an increase in 
As we will show later, both of these effects help to supress 
the bulk conductivity and stabilize the TI phase. 

The InAs/GaSb heterostructure we used comprises 
InAs and GaSb layers with nominal thicknesses of 12 and 
10 nm, respectively, sandwiched between Alo. 7 Gao. 3 Sb 
barriers [Fig. 2(a)]. The structure was grown by molecu¬ 
lar beam epitaxy on an n + -GaAs substrate which serves 
as a back gate. To bring the Fermi level close to the 



FIG. 2. (Color online) (a) Schematic of sample structure, (b) 
Six-terminal device pattern and typical nonlocal resistance 
measurement setup using dual lock-in technique. The dis¬ 
tance between adjacent contacts is 2 /im. The side contacts 
on the opposite edges (2-6 and 3-5) are separated by 2 /im. 


energy gap region, the upper Alo. 7 Gao. 3 Sb barrier is <5- 
doped with [Be] = 5 x 10 11 cm -2 at a setback of 5 nm .- 24 
A small six-terminal device with AuGe Ohmic contacts 
[Fig. 2(b)] is formed using deep ultraviolet lithography 
and wet etching. The distance between adjacent contacts 
is 2 /im. A Ti/Au front gate is evaporated on a 20-nm- 
thick AI 2 O 3 gate insulating layer deposited by atomic 
layer deposition^ 

Transport measurements were performed using a lock- 
in technique at a temperature of T = 0.27 K, unless 
otherwise noted. For nonlocal resistance measurements, 
alternating currents with different frequencies (19 and 23 
Hz) are injected from different contacts [Fig. 2(b)]. For 
each pair of voltage probes, two lock-in amplifiers are 
used to monitor the two frequency components. This 
dual lock-in setup allows us to simultaneously measure 
nonlocal resistances for different current injection paths. 
As we will see later, this method is particularly useful for 
the present purpose: it not only saves measurement time, 
but also allows for rigorous comparison between results 
for different current paths without its being influenced by 
irreproduciblc time-dependent fluctuation of the sample 
characteristics. 

Figure 3(a) shows longitudinal resistance ^ 14,23 at dif¬ 
ferent back-gate voltages (Vbg) plotted as a function 
of front-gate voltage (Vfg)- The suffixes of Rij.kt (= 
Vki /Iij ) indicate the contacts used for driving the current 
(*, j) and measuring the resultant voltage ( k,l ). Because 
of the Be doping in our sample, the Fermi level is in the 
valence band at Vfg = 0 V. With increasing Vfg, the 
Fermi level moves toward the conduction band. The re¬ 
sistance curve shows a peak when the Fermi level passes 
through the gap region. When the back gate is unbiased 
(Vbg = 0 V), the peak resistance is low (< 1 kfl). As 
Vbg is increased, the resistance peak shifts to lower Vfg 
and grows noticeably, reaching ~ 30 kfl at Vbg = 10 V. 
At the same time, spike-like fluctuations superimposed 
on the peak become conspicuous. These fluctuations are 
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FIG. 3. (Color online) (a) Longitudinal resistance (I?i4,23) 
at different back gate voltages (Vbg) as a function of front 
gate voltage (Vfg). Inset: magnified view of the the J?i 4,23 
traces in the low-resistance regime, (b) Peak value of f?i 4 : 23 
obtained from Gaussian fitting of the f?i 4,23 traces in (a), 
plotted as a function of Vbg- (c) Temperature dependence 
of -Ri 4,23 peak at different Vbg- The data are plotted as a 
function of inverse temperature (1/T). The short-dashed line 
is the Arrhenius fit to the data at Vbg = 10 V, with the slope 
corresponding to a gap of 0.2 meV. 
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mostly reproducible for repeated gate sweeps, except for 
irreproducible slow time-dependent onesJ^ 

The height of the # 14,23 peak deduced by fitting it 
with a Gaussian is plotted in Fig. 3(b) as a function of 
Vbg<^ The horizontal dashed line in the figure indicates 
the value, h/ 2e 2 (~ 13 kfl), expected from the quan¬ 
tized edge transport for the geometry used (h: Planck’s 
constant, e: elementary charge). At low Vbg> the 
peak resistance is much lower than this value, indicat¬ 
ing that considerable bulk conduction exists even when 
the Fermi level is in the gap region. As Vbg is in¬ 
creased, the peak resistance increases significantly, ex¬ 
ceeding h/ 2e 2 at Vbg ~ 7 V, and then tends to satu¬ 
rate at higher Vbg- Although edge transport in a 2 D 
TI should be protected from backscattering, in actual 
systems the conductance quantization is disrupted by 
processes that equilibrate the counterpropagating edge 
channels, as seen from the length dependence of the edge 
conductance^^ The most probable cause is the exis¬ 
tence of electron and hole puddles coming from the spa¬ 
tial potential fluctuations ) 27 ! 28 which can terminate the 
edge transport even in small samples and make the resis¬ 
tance higher than the expected quantized value d 0 ! 13 ! 29 ! 30 


FIG. 4. (Color online) (a)-(d) Nonlocal resistances (a) 
Ai6,34, (b) Ai6,45, (c) #21,34, and (d) #21,45, measured at 
Vbg = 0 and 10 V as a function of Vfg- (e) and (f) Ra¬ 
tios of nonlocal resistances (voltages) measured with adjacent 
voltage-probe pairs. Results for different current paths are 
compared for (e) Vbg = 10 and (f) 0 V. (g) and (h) Equiva¬ 
lent circuits for TI and normal conductor, respectively. 


Figure 3(c) shows the temperature dependence of the 
i?i 4_23 peak, plotted as a function of inverse temperature. 
At Vbg < 2 V, the peak resistance is almost temperature 
independent, indicating that the system does not have a 
full energy gap. At Vbg = 4 V the resistance is thermally 
activated in the high-T regime, suggesting the opening 
of an energy gap. However, the temperature dependence 
saturates at low-T regime (< 1 K), and the system re¬ 
mains fairly conductive, with the resistance well below 
h/ 2e 2 . With increasing Vbg further, the activated be¬ 
havior in the high-T regime becomes more pronounced^ 
and the resistance value in the low-T regime increases. 

In our previous study, we compared samples with dif¬ 
ferent InAs layer thicknesses (w = 10, 12 , and 14 nm) 
while keeping Vbg = 0 The temperature depen- 
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dence for Vbg = 2 and 6 V shown in Fig. 3(c) is similar, 
respectively, to that of the w = 12 and 14 nm samples (at 
Vbg = 0 V) in the previous report, which we identified as 
being a TI (w = 12 nm) and a semimetal ( w = 14 nm). 
As seen in Fig. 3(c), the resistance increases further at 
higher Vbg (> 8 V). However, the behavior at Vbg = 8 
and 10 V is distinct from that of the normal band insu¬ 
lator previously reported for the w = 10 nm sample (at 
Vbg = 0 V), in which the resistance increases to ~Mf2. 
Here, it must be noted that the sample studied here was 
fabricated from a wafer different from the one used for 
the w = 12 nm sample in the previous study. Compar¬ 
ing the results for the two samples with w = 12 nm, we 
infer that the actual InAs layer thickness is larger in the 
sample used in this study. 

The results in Fig. 3(c) show that the contribution of 
bulk conduction can be reduced by applied gate voltages. 
The different contributions of bulk and edge transport at 
low and high Vbg can be further highlighted by nonlo¬ 
cal measurements. Figures 4(a)-(d) display nonlocal re¬ 
sistances measured in four different configurations. The 
black and red curves in each panel represent the results 
obtained for Vbg = 0 and 10 V, respectively. The four 
traces for each Vbg were all measured simultaneously us¬ 
ing the dual lock-in technique. The horizontal dashed 
line in each panel indicate the value, hj 6 e , expected for 
the quantized edge transport in the geometry used for 
the nonlocal measurements^ The nonlocal resistances 
for Vbg = 0 V (which are enlarged in the figures by a 
factor of 20) are far below this value, h/ 6 e 2 . Similar 
to the case of longitudinal resistance, the nonlocal re¬ 
sistance peak grows significantly when Vbg = 10 V is 
applied, with the peak values exceeding hj 6 e 2 . 

As shown in Ref. [U that the transport is dominated 
by edge channels can be confirmed by examining if the 
ratio of the nonlocal resistances (or voltages) measured 
with adjacent voltage-probe pairs depends on the current 
path. The results for Vbg = 10 V are shown in Fig. 4(e), 
where we compare the voltage ratios V 45 /V 34 for two cur¬ 
rent paths Il6 and I 21 as a function of Vfg- In the Vfg 
range corresponding to the peak in the nonlocal resis¬ 
tance (0.3 < Vfg < 0.7 V), the voltage ratios for the two 
current paths completely agree with each other. Such 
agreement has been confirmed for all contact geometries 
[Fig. 5(d) and Fig. 6 (Appendix)]. As shown in Ref. [H, 
these results indicate that the current flows only along 
the edges of the sample, as illustrated in Fig. 4(g). This 
can be understood by noting that if the current also flows 
through the bulk region, the voltage distribution in the 
sample should depend on the current path [Fig. 4(h)]. In 
fact, the latter behavior is seen in the low- and high- Vfg 
regions of Fig. 4(e), where the voltage ratios no longer 
agree with each other. This shows that in these Vfg 
ranges the central bulk region becomes conductive, with 
the majority carrier type’s being p and n on the low- and 
high-VFG sides, respectively. We note that also in these 
Vfg ranges the voltage ratio depends on VfGj indicat¬ 
ing that edge transport coexists with bulk conduction. 
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FIG. 5. (Color online) Voltage ratios V 23 /V 34 (= 

Rij, 23 /Rij, 3 i) for current paths I\e and 765, measured at (a) 
Vbg = 0, (b) 4, (c) 6, (d) 10 V. 


Similar results were obtained from nonlocal resistance 
measurements at Vbg =8 V. 

At Vbg = 0 V [Fig. 4(f)], the voltage ratios are only a 
weak function of Vfg> and different current paths yield 
different voltage ratios. Similar results were confirmed 
for other contact geometries. We note that when the 
bulk transport is dominant the voltage distribution in 
the sample is determined solely by the contact geome¬ 
try and does not depend on sample conductivity. Thus, 
the results in Fig. 4(f) indicate that at Vbg = 0 V the 
transport is dominated by bulk conduction even when 
the Fermi level is in the expected gap region, implying a 
semimetallic band structure. 

We study the evolution from a semimetal to a TI in 
more detail using nonlocal measurements. Figures 5(a)- 
(d) display how the voltage ratios V 23 /V 34 for two current 
paths 7i6 and evolve with Vbg- With increasing Vbg , 
the voltage ratios acquire Vfg dependence, approaching 
each other, and at Vbg = 10 V they completely overlap in 
the range of 0.4 < Vfg 0.7 V. At Vbg = 4 V, the volt¬ 
age ratios do not agree, indicating a finite contribution of 
bulk conductivity. This is consistent with the longitudi¬ 
nal resistance peak’s being lower than h/ 2e 2 at Vbg = 4 
V [Fig. 3(c)]. The thermal activation of the longitudinal 
resistance in the high-T regime, in turn, indicates that 
at Vbg = 4 V the bulk conduction occurs locally. This 
implies that semimetallic and TI regions coexist within 
the sample when the Fermi level is in the gap region. 
Spatial variation of the carrier density or layer thickness 
can cause such sample inhomogeneity. 

Knez et al. also examined the effects of an electric field 
on residual bulk conductivity gtmik in GaSb/InAs de- 






5 


vices and observed that gtmik decreased with applied gate 
voltagesi 12 ' 32 However, in their undoped samples (/bulk re¬ 
mained significant (> 2 e 2 /h) even when the largest elec¬ 
tric field was applied^ The results were interpretted in 
terms of the theory in Ref. [H, which states that even 
weak scattering can destroy the anticrossing gap and re¬ 
sult in finite bulk conductivity. For T <C A <C -Eho — E e o 
(where T is the level broadening), the theory provides 
a crude estimate of the residual conductivity: <?buik ~ 
( e 2 /h)(E h0 — E e0 )/A. The effect of the electric field is 
then understood to decrease (E h0 — E e0 )/ A. However, 
the theory does not account for the vanishing of bulk con¬ 
ductivity we observed at high Vbg- This suggests that 
the assumption T<A<C -Eho — E e0 no longer holds in 
the high-V bg regime of our experiments, where E\^ — E b q 
could be comparable to A depending on the InAs layer 
thickness. 

Another possible mechanism that gives rise to resid¬ 
ual bulk conductivity is the anisotropy of the GaSb va¬ 
lence band^ In the presence of band anisotropy, the 
energy position of the anticrossing gap depends on the 
k direction. If the anisotropy is large as compared to 
A, the anticrossing gaps for different k directions may 
not overlap in energy, resulting in a semimetallic band 
structure [Figs. 1(c) and (e)]. The effects of the electric 
field are to shift the anticrossing to a smaller k, where 
the effect of the GaSb valence band anisotropy is smaller 
[Fig. 1(d) and (f)], and enlarge A/ 7 Both effects would 
work to suppress the bulk conductivity and stabilize the 
TI phase. Most importantly, our nonlocal measurements 
clearly show that the bulk conduction vanishes at higher 
Vbg (= 8 and 10 V), demonstrating that an originally 
semimetallic heterostructure can be tuned into a TI. Al¬ 
though further studies are necessary to clarify the role of 
disorder; 16 ! 34 our results clearly represent an important 
aspect of the InAs/GaSb system, which will be useful for 
the study of TIs. 

In summary, we have shown that a semimetallic 
InAs/GaSb heterostructure can be tuned into a TI by ap¬ 
plying gate voltages. Nonlocal resistance measurements 
using a dual lock-in technique enabled us to confirm 
the suppression of bulk conductance and the transition 
to a topological insulator without influence from time- 
dependent resistance fluctuations. Our results warrant 
further investigations to explore the transition between 
nontopological and topological insulators. 


ACKNOWLEDGMENTS 


We sincerely thank Y. Ishikawa and H. Murofushi for 
sample preparation. This work was supported by JSPS 
KAKENHI Grant Number 26287068. 


Appendix 

Figure 6 summarizes the results of nonlocal resistance 
measurements at Vbg = 10 V for contact geometries 
other than those in Figs. 4(e) and 5(d). For all the geome¬ 
tries, the nonlocal resistance ratios for different current 
paths precisely match each other in the shaded region, 
demonstrating that the bulk region is insulating and the 
transport is governed by edge channels in the correspond¬ 
ing Vfg range. 
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FIG. 6. (Color online) Nonlocal resistance measurements 
at Vbg = 10 V for various contact geometries shown in the 
insets. In each panel, nonlocal resistance ratios for different 
current paths are plotted as a function of Vfg (T = 0.27 K). 
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